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Abstract 

In the present paper, the connection between surface order-disorder phase transitions 
and the percolating properties of the adsorbed phase has been studied. For this purpose, 
four lattice-gas models in presence of repulsive interactions have been considered. Namely, 
monomers on honeycomb, square and triangular lattices, and dimers (particles occupying 
two adjacent adsorption sites) on square substrates. By using Monte Carlo simulation and 
finite-size scaling analysis, we obtain the percolation threshold 9c of the adlayer, which 
presents an interesting dependence with w/fc^T (being ks and T, the lateral interaction 
energy, the Boltzmann's constant and temperature, respectively). For each geometry and 
adsorbate size, a phase diagram separating a percolating and a non-percolating region is 
determined. 
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I. INTRODUCTION 



Despite over three decades of intensive work, the interplay between percolation prop- 
erties and thermal phase transitions is still an open problem. In this sense, the study 
of geometrical structures close to the critical point allows a better understanding of the 
mechanism of the phase transitioo^'^'^. The geometric critical phenomena exemplified by 
percolation possess many striking parallels with the thermally driven critical phenomena, 
as it is provided by the Fortuin-Kasteleyn mapping^. Cluster description of thermody- 
namic phase transitions have been used since long time to elucidate the nature of the 
transitions by providing a geometrical interpretation of density correlations^. Fisher has 
introduced the phenomenological droplet model in which the fluctuations are associated to 
clusters or droplets which percolate at critical point^. More recently, it was concluded that 
the percolation transition can be considered as a particular case of the g-state Potts model 
where q is equal to l^'^, and can be described as a second order one, whose universality 
class depends only upon the space dimensionality. By following this line of reasoning, a 
wide variety of systems have been studied. Among the more recent contributions, the 
behavior of colloids and gels has been discussed in an interesting paper by A. Coniglio^. 
The author concluded that it is very important to define the appropriate cluster for each 
phenomenon. Later, an important contribution has been made by S. Fortunato^*^ who 
analyzes the critical exponents of both the thermal and percolation phase transitions 
occurring for different models in two dimensions. 

Although most of the works on the subject are devoted to the study of lattice-gas models 
in presence of attractive lateral interactions (or ferromagnetic coupling in magnetic lan- 
guage), there have been a few studies related to repulsive interactions and order-disorder 
phase transitions. In the last case, the definitions of connectivity and cluster of particles 
belonging to the adsorbed phase are not the same and the relationship between perco- 
lating clusters and critical points is non-trivial. In a previous paper—, we studied the 
percolation of monomers on a square lattice as the particles interact with repulsive ener- 
gies. The present contribution goes a step further, including honeycomb and triangular 
substrates and multisite-occupancy^'i^i'i^'i^'i^ (adsorbates occupying more than one site). 

The outline of the paper is as follows: In Section [TTl we describe the lattice-gas model 
and the finite-size scaling theory. In Section IIIII we present and discuss the results along 



with general conclusions. 



II. MODEL AND FINITE-SIZE SCALING THEORY 

In order to consolidate the ideas involved here, four different physical systems have been 
considered, according to the adsorbate's size and surface geometry: 
Model I): Monomers adsorbed on square lattices. 
Model II): Monomers adsorbed on honeycomb lattices. 
Model III): Monomers adsorbed on triangular lattices. 
Model IV): Dimers adsorbed on square lattices. 

In all cases, the substrate is represented by M = L x L equivalent adsorption sites with 
periodic boundary conditions. 

To describe the system of particles adsorbed on M sites at a given temperature T, 
let us introduce the occupation variable q which can take the following values: q = if 
the corresponding site is empty and q = 1 if the site is occupied by an adatom (or dimer 
unit). Under these considerations, the Hamiltonian of the system is given by, 

M 

H = wJ2ciCj- N{k-l)w + eoJ2ci (1) 

(i.j) i 

where w is the nearest-neighbor (NN) interaction energy (we focus on the case of repulsive 
lateral interactions among adsorbed particles, w > 0); represents pairs of A^A^ sites 
and k = 1(2) for monomers(dimers). The term N{k — l)w is subtracted in eq. since, 
in the case of k = 2, the summation over all the pairs of NN sites overestimates the 
total energy by including bonds belonging to the A^ adsorbed dimers. Finally, eo is 
the adsorption energy of the sites on the surface (we have taken = without loss of 
generality) . 

For fixed values of surface coverage, 6 = kN/M, and temperature T, the thermody- 
namic equilibrium is reached in the canonical ensemble by using a standard Kawasaki 
algorithn>iiii^. Thus, a set of m = 2000 samples in thermal equilibrium is generated by 
taking configurations separated from each other by 1000 Monte Carlo steps in order to 
avoid memory effects. 

The central idea of the percolation theory is based on finding the minimum concentration 
6 for which at least a cluster [a group of occupied sites in such a way that each site has at 
least one occupied nearest-neighbor site] extends from one side to the opposite one of the 



system. This particular value of the concentration rate is named critical concentration or 
percolation threshold and determines a phase transition in the system. 

It is well known that it is a quite difficult matter to analytically determine the value of 
the percolation threshold for a given lattice^'^'^'^'^. Thus, in most cases, percolation 
thresholds have to be estimated numerically by means of computer simulations. 

As the scaling theory predicts^, the larger the system size to study, the more accurate 
the values of the threshold obtained therefrom. Thus, the finite-size scaling theory give us 
the basis to achieve the percolation threshold and the critical exponents of a system with a 
reasonable accuracy. For this purpose, the probability R = Ri{9) that a lattice composed 
of LxL elements (sites or bonds) percolates at concentration 6 can be defined^. Here, the 
following definitions can be given according to the meaning ofX: a) rI^^\9) = the prob- 
ability of finding a rightward (downward) percolating cluster; b) Rl{0) = the probability 
that we find a cluster which percolates both in a rightward and in a downward direction; 
c) R^{0) = the probability of finding either a rightward or a downward percolating clus- 



ter and d) Rt{9) = i [r^{9) + RE{0)] = \ {R^A&I 
and using the methodology described in Refs. U 
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colation thresholds were calculated by extensive use of finite size scaling techniques. The 
interested reader is urged to read the above cited articles for a more complete discussion 
of this issue. 



III. PERCOLATION PHASE DIAGRAM: RESULTS AND CONCLUSIONS 

By using the scheme discussed above, the critical curves, Qc vs. K (being K = w/kBT), 
separating the percolating and non-percolating regions, were calculated. 

In Fig. 1, the percolation phase diagram is shown for Model I. As = (non-interacting 
adsorbate), the adsorption-desorption process reproduces a random deposition, which is 
fully equivalent to random percolation. Consequently, the expected value of 6c = 0.592 is 
reached for = 0. As -ft" is increased, two well-differentiated regimes can be distinguished: 
i) from i^T = up to i^' ~ 1.76 (being K ^ 1.76 the reduced critical temperature for 
the order-disorder phase transition occurring in the system), 6c increases linearly with 
K; and ii) for K > 1.76, 6c remains constant as K is increased. This behavior can 
be explained from simple geometrical arguments. Namely, lateral repulsive couplings 
avoid the occupation of nearest-neighbors sites, and consequently, increase the percolation 



threshold. In the hmit case, once Kc is reached, the adlayer does not vary significantly 
as K is increased, and 9c reaches its saturation value, being O^. 0.66 for K > 1.76. It 
is worth to emphasize that the presence of strong lateral interactions (and consequently, 
the existence of a phase transition occurring in the system) yields an increase in the 
computational effort to get accuracy values of the percolation threshold. 

It is important to bear in mind that the points in Fig. 1 correspond to states in thermal 
equilibrium. In order to reflect this situation, we have calculated the adsorption isotherms 
(mean coverage as a function of the reduced chemical potential, /i/fc^T) for repulsively 
interacting adparticles in a wide range of temperatures. The adsorption process was simu- 
lated through a Grand Canonical Ensemble Monte Carlo (GCEMC) method. Relaxation 
toward equilibrium relied upon Glauber dynamics^-. 

For a given value of temperature T and chemical potential /i, an initial conflguration 
with N monomers (dimers) adsorbed at random positions on (2A^) sites is generated. 
Then an adsorption-desorption process is started, where a site (pair of nearest-neighbor 
sites) is chosen at random and an attempt is made to change its occupancy state with 
probability given by the Metropolis rule^: 

P = min|l,exp 

where H = Hj — Hi is the difference between the effective Hamiltonians of the flnal and 
initial states, being H = H — fiJ^Ci- A Monte Carlo Step (MCS) is achieved when 
sites (pair of sites) have been tested to change its occupancy state. The equilibrium state 
can be well reproduced after discarding the flrst m' = 10^ — 10^ MCS. Then, averages are 
taken over m = 10^ — 10^ successive configurations. In this framework, the mean coverage 
is obtained as: 

1 ^ 

^ = ^E<c.> (3) 

i 

where the thermal average (...), means the time average over the Monte Carlo simulation 
run. 

In order to compare our numerical results with a theoretical prediction, we have used 
one of the most reliable methods for studying the thermodynamic properties of a system 
suffering a phase transition: the Real Space Renormalization Group (RSRG)^. The 



interested reader is referred to Ref. 
and to Refs. 38 39 3: 



36 



37| for a detailed description of the RSRG method 



41| for applications of the RSRG method to lattice gas models. 



In the RSRG method developed by Niemeyer and van Leeuwen^ and Nauenberg and 
Nienhuis^'^, the whole lattice is divided into blocks (or cells) of L sites. A block spin iSq, 
is assigned to each block. All blocks together must form a square lattice with the lattice 
constant \/La. The RSRG transformation of the spin system allows the reduction of the 
number of independent variables, i.e. the transition from the set of N site spins {sj} to 
N/L block spins {S^j-We note that two values of the block spin = ±1 corresponds to 
2^ site spin configurations (since L spins are combined to form a block). For blocks with 
odd number of spins is usually determined by the so-called "majority rule" (MR)^. 
For even L a rule must be introduced in order to assign a definite value of the block 
spin to any given configuration with the sum of site spins equal to zero. In any case an 
obvious condition must be fulfilled: if the site spin configuration {si, S2 ■ ■ ■ sl] is assigned 
to a block spin with weighting factor P, then the configuration {— Si, —S2 . . . — si} is 
assigned to the —S^ domain with the same P. 

In the framework of the RSRG approach, one usually employs periodic boundary con- 
ditions. It is assumed that the whole lattice is given by the periodic continuation of a 
small cluster of blocks. In our calculations we consider the smallest possible cluster of 
two blocks. Due to the simplicity of this cluster, no additional interactions appear in 
the renormalized Hamiltonian. It is the same Hamiltonian of the square Ising spin sys- 
tem with, however, renormalized values for the external magnetic field and for the pair 
interaction parameter. 

As was shown by Nauenberg and Nienhuis^, the free energy of the system for any values 
of magnetic field and interaction parameter can be evaluated in a series of sequential RSRG 
transformations of the original Hamiltonian. The interested reader will find details of the 
application of this technique for Models I, II and III in the following references: square 



lattice, Ref. 



45], triangular lattice, Ref. |46] and honeycomb lattice Ref. 47]. 



The results, obtained by Real Space Renormahzation Group, RSRG (solid lines) and 
Monte Carlo, MC, methods (small squares), are shown in Fig. 2^-. At high temperatures 

the isotherms are close to the Langmuir case (lattice-gas without lateral interaction), i.e. 

exp/3(/x + £o) 
'^'^ = l + exp/5(/. + ..)- 
At low temperatures a broad plateau occurs around half coverage. This plateau corre- 
sponds to the c(2 X 2) ordered lattice-gas phase (or, in magnetic language, to the AF 
ordered two-dimensional spin system). Large spheres in Fig. 2 are the same points as in 



Fig. 1. It is clear that such hne avoids to enter in the region of the coexistence of phases. 

To reinforce the above result, in Fig. 3 the percolation line is plotted (spheres) together 
with the coexistence curve in the temperature-concentration diagram, which limits the 
region where the c(2 x 2) ordered phase percolates. As it can be observed, the percolation 
line remains in the region where the system is disordered. These features clearly reveal 
that the definitions of connectivity (in the sense of standard random percolation) and 
therefore the definitions of the clusters of c(2 x 2) ordered structure are not the same. 

In the case of Model II (Fig. 4), the general trend is similar to that of the square 
lattice. The curve grows monotonically up to a value of ~ 2, where it reaches an 
almost constant value of 0.75 for the critical coverage degree. The explanation of this 
trends is similar to the first case. The percolation line (spheres) plotted together with 
the adsorption isotherms for repulsively interacting particles adsorbed on a honeycomb 
lattice is shown in Fig. 5. In Fig. 6 the same line is plotted in conjunction with the phase 
diagram. Again, the percolation line remains in the region where the system is disordered. 

For Model III (Fig. 7), it can be observed that the value of the percolation threshold is 
near 0.5 in the whole range of K (notice the scale in the graph). A complete understand- 
ing of the phase diagram is a very important help in the description of the peculiarities of 
the temperature dependence of the percolation threshold. In order to explain the antifer- 
romagnetic ordering we recall that a triangular lattice can be seen as a system composed 
of three equivalent triangular sublattices. As is well known, pairwise interaction results 
in a symmetrical phase diagram around 9 = 0.5. For this lattice-gas system, triangular 
antiferromagnetic lattice-gas, the phase diagram (T, 6) consists of two symmetrical curves 
around 6 = 0.5. For 6 < 1/3, the ordered phase reveals that the particles are arranged in 
such a way that pairs of particles on nearest-neighbor lattice sites are not present. In fact, 
most of the adsorbed particles are located in only one of the three equivalent triangular 
sublattices, thus avoiding possible interactions with other particles (which is equivalent 
to an ordered t|| phase using magnetic language). This ordered phase prevails over the 
range 6' < 0.5 as is indicated in Fig. 8. 

The symmetric branch of the phase diagram reflects the ordered phase where two of the 
sublattices are occupied (HJ, in the magnetic language). This phase diagram has been 
investigated by Schick, Walker and Wortis^*^ by using RSRG as well as by the transfer 
matrix method^ and MC simulations^*^. 



It is interesting to note that for 9 < 0.5 a) the corresponding ordered phase percolates 
but b) there is not standard percolation of the adsorbed monomers on the lattice. The 
percolation line dividing the percolating and non percolating area yields in the disordered 
phase of the antiferromagnetic phase diagram (see Fig. 8). Furthermore, in the whole 
range of temperature, the percolation occurs for a coverage lower than what is needed to 
built the ordered phase for ^ > 0.5 

Fig. 9 shows the percolation line together with the adsorption isotherms for repulsively 
interacting particles. The adsorption isoterms present clearly defined plateaus located at 
= \ and = §• Strong enough repulsion produces ordered phases when particles occupy 
preferentially sites of a single sublattice (6* = |), or two sublattices {6 = |). 

In the case of dimers on square lattices (Fig. 10), the curve 9^ vs. K (spheres) is 
similar to that of Fig. 1 (open circles). Thus, 9c grows linearly from 9c{K = 0) = 
0.56 (as it is expected for random percolation of dimers) to a saturation value close 
to 2/3. This behavior has interesting consequences on the temperature-concentration 
phase diagram. In fact, as it has been reported in the literature^, a "zig-zag" (ZZ) 
ordered phase, characterized by domains of parallel ZZ strips oriented at ±45° from 
the lattice symmetry axes, separated from each other by strips of single empty sites, 
was found at 2/3 monolayer coverage (Fig. 11). The ordered phase is separated from 
the disordered state by a order-disorder phase transition occurring at a finite critical 
temperature. An accurate determination of this critical temperature has been recently 
obtained [kBTjw = 0.182(1)]^. 

A simple inspection of Fig. 11 shows the existence of long-range connectivity for the 
low-temperature phase at 2/3 coverage. This finding, along with the tendency to 2/3 of 
the curve in Fig. 10, clearly reveals the interplay between the surface order-disorder phase 
transition and the percolating properties of the adsorbed phase at 2/3 monolayer cover- 
age. Namely, i) the ZZ ordered phase represents the state of the adlayer at percolation 
threshold and K ^ oo and ii) the curve 9c vs. K crosses the coexistence line on the 
temperature-concentration phase diagram at {9c = 2/3, fc^Tc/w = 0.181) and penetrates 
in the ZZ-region. The last point represents the main difference between Model IV and 
the other models previously analyzed. A systematic analysis of critical exponents was not 
carried out since this was out of the scope of the present work. 

Since properties of adsorbed layers are often determined by measuring adatom coverage 



versus adsorbate gas pressure, it is important to show what adsorption isotherms would 
look in this case. Thus, Fig. 12 shows a set of adsorption isotherms for dimers and 
different values of repulsive nearest-neighbor interactions together with the percolation 
line. Since the symmetry particle-vacancy, valid for monoatomic species, is broken for 
dimers, the adsorption isotherms are not symmetric around ^ = 0.5. In addition, two well 
defined and pronounced steps appears as K increases. At 9 = 0.5, a well defined array of 
dimers resembling a c(2 x 2) phase, is found. The ordered structure is characterized by 
a repetition of alternating files of adsorbed dimers separated by 2 adjacent empty sites. 
As the chemical potential /i increases and 9 approaches 9 = 2/3, incoming dimers are 
adsorbed forming domains of parallel zig-zag rows (ZZ phase) as it was discussed above. 
These structures are clearly evidence of a low temperature ordered phase. In fact, the 
systems undergoes continuous phase transitions, from disorder to ordered structuresi^'^^'^. 

In summary, we presented a model to investigate the process of adsorption of interact- 
ing monomers on square, honeycomb and triangular lattices and studied the percolating 
properties of the adsorbed phase. By using Monte Carlo simulation and finite-size scal- 
ing theory, we obtained the percolation thresholds for different values of concentration 
and temperature. From this analysis, a critical curve in the 9 — T space was addressed. 
The line separating the percolating and non-percolating regions was explained in terms 
of simple considerations related to the interactions present in the problem. 
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Figure Captions 



Fig. 1: Phase diagram, 9^ versus which shows the curve separating the percolating 
and nonpercolating regions for the case of monomers on a square lattice. The vertical 
dashed line at X = 1.76 denotes the reduced critical temperature for the phase transition 
occurring in the adlayer phase for repulsive interacting particles. Horizontal dashed line 
at 9(. = 0.662 is the critical coverage at saturation regime for K > 1.76. The error bars 
are smaller than the symbol size. 

Fig. 2: Adsorption isotherms for Model I (surface coverage, 9, vs. normalized chemical 
potential, ^/ksT and reciprocal temperature expressed in units of K). Solid lines are 
obtained by the RSRG method, small symbols denote MC data while large spheres are 
the same points as in Fig. 1. 

Fig. 3: Phase diagram (critical temperature versus surface coverage) corresponding to 
Model I obtained by RSRG, solid line. Phase diagram, versus 9ci which shows 

the curve separating the percolating and nonpercolating regions, spheres. The inset is a 
snapshot of the ordered phase c(2x2). 

Fig. 4: Phase diagram, 9c versus K, for Model II. The critical coverage at saturation 
regime is 9c — 0.758. The error bars are smaller than the symbol size. 
Fig. 5: Adsorption isotherms for Model II and different values of K, as indicated with 
small symbols. Solid lines are obtained by the RSRG method, small symbols denote MC 
data while large spheres are the same points as in Fig. 4. 

Fig. 6: Phase diagram (critical temperature versus surface coverage) corresponding to 
Model II obtained by RSRG, solid hue. Phase diagram, K"^ versus 9c, which shows 
the curve separating the percolating and nonpercolating regions, spheres. The inset is a 
snapshot of the ordered phase. 

Fig. 7: Phase diagram, 9c versus K, for Model III. The error bars are included in the 
figure. 

Fig. 8: Phase diagram (critical temperature versus surface coverage) corresponding to 
Model III obtained by RSRG, solid line. It consists of two symmetrical curves around 
9 — 0.5. The first one, for 9 < 0.5 limits the ordered phase where the particles are 
arranged in such a way that pairs of particles on nearest-neighbor lattice sites are not 
present (which is equivalent to an ordered ][[ phase using magnetic language). The 



second one, reflects the ordered phase where two of the sublattices are occupied (tti)- 
The insets are snapshots of the corresponding ordered phases. Phase diagram, versus 
9c, which shows the curve separating the percolating and nonpercolating regions, spheres. 
Fig. 9: Adsorption isotherms for Model III and different values of K, as indicated with 
small symbols. Solid lines arc obtained by the RSRG method, small symbols denote MC 
data while large spheres are the same points as in Fig. 7. 

Fig. 10: Phase diagram 9^ vs. K for Model IV (spheres) plotted together with data of 
Model I for comparison (open circles). 

Fig. 11: Phase diagram (critical temperature versus surface coverage) corresponding to 
Model IV, solid line. It consists of curves representing the regions where the phases plotted 
in the insets are stables. Phase diagram, K^^ versus 9c, which shows the curve separating 
the percolating and nonpercolating regions, spheres. 

Fig. 12: Adsorption isotherms for Model IV and different values of K, as indicated. Small 
symbols denote MC data while large spheres are the same points as in Fig. 10. 
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